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Abstract 



It has been suggested by Roukema and coworkers (hereafter R04) that the topology of the Universe as probed by the "matched 
circles" method using the first year release of the WMAP CMB data, might be that of the Poincare dodecahedral space (PDS) model. 
An excess in the correlation of the "identified circles" was reported by R04, for circles of angular radius of ~ 11° for a relative phase 
twist —36°, hinting that this could be due to a Clifford translation, if the hypothesized model were true. R04 did not however specify 
the statistical significance of the correlation signal. 

We investigate the statistical significance of the signal using Monte Carlo CMB simulations in a simply connected Universe, and 
present an updated analysis using the three-year WMAP data. We find that our analyses of the first and three year WMAP data 
provide results that are consistent with the simply connected space at a confidence level as low as 68%. 

Key words, cosmology: observations - cosmology: cosmic microwave background - cosmology: cosmological parameters 



1. Introduction 

If the topology of the Universe were multiply connected, as op- 
posed to simply-connected, and if the comoving size of the fun- 
damental domain (FD) were smaller than the comoving distance 
to the surface-of-last-scattering (SLS), then it should be pos- 
sible to detect repeating patterns in the CMB fluctuations us- 
ing full-sky data of sufficient signal-to-noise ratio. These fluc- 
tuations would be those lying along pairs of circles defined by 
points of interse ction between differen t copies of the SLS in the 
covering space dCornish et alJll99"8bh . These patterns although 
found in different directions of the sky, would constitute so- 
called "matched circles", as they would represent the same phys- 
ical points, but observed from different directions due to topo- 
logical lensing. 

While this principle is true for any 3-manifold model of space, 
the number of pairs of "matched circles" or their sizes and rela- 
tive spatial orientations, as well as their handedness, or phase 
shift, depend significantly on the assumed 3-manifold and its 
topological properties, thus providing a way to observationally 
distinguish between models. 

While the positive correlation signal from matched pairs is 
expected dire ctly from the metric p erturbations, via the Sachs- 
Wolfe Effect dSachs & Wolfelll967h . there are many other cos- 
mological effects (e.g. t he Doppler effect, the Integrated Sachs- 
Wolfe E ffect (ISW)) (iKunz et alj |2008|). astrophysical fore- 
grounds (iBennett et al.ll2003 ) and instrumental effects that con- 
stitute noise, from the point-of-view of a matched circles search, 
and the magnitude of the effects depends on the angular scale. 

Although the CMB data have been analyzed to de- 
tect topologi cal lensing sigii als since the availability of the 
COBE data (! Roukem a| l20Ciol) . the release of the WMAP ob- 
servations has provided full-sky data of unprecedented ac- 
curacy and resolution, opening up more promise for di- 
rect tests of the topology of the Universe. Although the 



"matched circles" test is straightforward, it is limited due to 
noise and FD size constraints. Additional theoretical predic- 
tions can be used as independent tests that involve predic- 
tions of CMB temperature and polarization fluctuations for 
the case that the Universe is multiply connected, both in 
real and spherical harmonic spaces, or topologi cal effects on 
the CMB power spectruni dAurich et al. 2004 IWeeksl 120031: 



Gomero & Reboucas '2003!; pineen et al.l i2 005: Riazue lo et al 
2004t iLehoucq et al.j jl999ninouel 119991; ,Niarchou & Jaffe 



2OO7I: ICornish et al.l H998al; Ide Oliveira-Costa & SmootI 1 1995 



Kunz et al."2006'; 'Luminet et al."2003'; 'Niarch ou & Jaffel bood; 
Phillips & Kogut 2006; Caillerie et al. 2007). Although a suc- 
cessful "matched circles" test would provide strong support 
for the Universe being multiply conne cted, no statis t ically- 
significant evid ence has been found dCornish et all 120041: 
iKev et al.ll2007h. 

In iRoukema et al.l d2004l) we performed a "matched 

circles" searc h using the first year WMAP ILC map 

dBennett et al.l l2003h and found an excess correlation, which 
one would expect under the PDS hypothesis for circles of 
angular radii a ~ 11° with centers towards {l°,b°) — 
{(252, 65), (51, 51), (144, 38), (207, 10), (271, 3), (332, 25)} 
and their opposites (Fig.[T]i. 

In this present work we have two key objectives. We revisit 
those results, verify the existence of the excess correlations and 
quantify their statistical significance. Secondly, we update the 
search with the WMAP three year data release, extend it to probe 
three different resolutions (smoothing lengths) and define the de- 
tection confidence thresholds. We also discuss the effects of un- 
derlying 2-point correlations, smoothing length and incomplete 
sky coverage on the value of correlation coefficient. 

In section [2] we introduce the datasets used in the analysis, 
provide details of their preprocessing, and describe simulations 
that we use to complete a statistical significance analysis. In sec- 
tion [3] we introduce the details of the statistics being performed 
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Fig 1. Visualization o f the matched circles solution reported in 
iRoukema et al.l (l2004h and reproduced in this paper to constrain 
its statistical significance, over-plotted on the first year ILC map 
masked with Kp2 sky mask. 

and our confidence-level analysis. Results are presented in sec- 
tionH] We conclude in section|6] 

2. Data and simulations 

We perform a "matched circles" search using two sets of data. 
Fi rstly, for the sake of co mpatibility, we choose the same data as 
in IRoukema et al.l ( |2004 - i.e. the first year WMAP ILC map. 

The topologically interesting signal generally dominates 
over the Doppler (a nd other) components on large scales 
jRiazuelo et al.l l2004b : this is a motivation for using a large 
smoothing length. However, extended flat fluctuations that hap- 
pen to have a similar lar ge scale trend can lead to false pos- 
itives on large scales (KexeLaL 2007). This implies a trade 
off in the choice of the smoothing length of the data, between 
large smoothing lengths preferred by the topologically interest- 
ing content, and small smoothing scales which avoid false pos- 
itives induced by chance correlations of extended flat fluctua- 
tions. 

We choose to test three different smoothing scales: 
FWHM = A e {1°,2°,4°}. 

The ILC map was obtained from a linear combination of one 
degree smoothed maps in the five frequency bands, by inverse 
noise co-adding them, hence its resolution is consistent with a 
one degree smoothing scale. We further Gaussian smooth this 
map in spherical harmonic space by convolving it with Gaussian 
beam response kernels of FWHM corresponding to 2° and 4° 
respectively, to obtain the first set of data for matched circles 
tests. 

Secondly, we choose the three year foreground reduced 
WMAP data from individual frequency bands Q[l/2],V[l/2] and 
W[ 1/2/3/4] and co-add them into one map, accordiiig to th e in- 
verse noise weighting scheme used in lHinshaw et al.N2003b . We 
call the resultant map the "INC map". We smooth the INC map 
using a Gaussian convolution kernel to four different FWHM 
smoothing lengths: A £ {0.5°, 1°,2°,4°}. 

We downgrade all data from t he initial resolution defined by 
the Healpix pixelization scheme dGorski et al.ll2005T) with reso- 
lution parameter Us = 512 (res. 9) to a resolution parameter of 
Us = 256 (res. 8). We remove the the residual monopole and 
dipole components {£ = 0, 1) in spherical harmonic spaceQ be- 
cause these components are of no cosmological interest. At the 

' Since the residual WMAP maps foregrounds are strong, we perform 
this step using the Kp2 sky mask to keep the compatibility between the 
data and simulations. 



final stage of preprocessing, we remove the residual monopole 
by offsetting the maps in real (2-sphere) space so that (T) = 
outside the Kp2 sky mask. 

Throughout the analysis, i.e. for both the ILC and INC maps, 
we use the Kp2 sky mask, which masks ~ 15% of the sky in- 
cluding the brightest resolved point sou rces. The Kp2 sky mas k 
is different from the sky mask used in IRoukema et al.1 (12004 . 
While the ILC map is best suited e.g. for the full sky low multi- 
poles alignment analysis, for the purpose of the matched circles 
test, the residual galactic contamination should be masked out, 
although we realize that the use of the Kp2 mask may be too 
conservative. In App. IA.2I we compare the impact of different 
sky masks and demonstrate that our results are not very sensitive 
to the precise characteristics of the sky mask. 

For each of the two data sets, we produce N^i^a = 100 realis- 
tic Gaussian random field (GRF) signal and noise simulations of 
the WMAP data to quantify the statistical significance of plausi- 
ble detections, to discard false positives, and to resolve the 2a- 
confidence levels. Therefore, for the first dataset we simulate the 
first year ILC map, inside "region 0" defined outside the Kp2 
sky mask of Bennett et al. (2003), and for the second dataset we 
simulate the three year WMAP INC map. 

As will be shown in Sect. |3] the matched circles correlation 
coefficient depends on the monopole value in the map. Also, 
in principle it is sensitive to the shape of the 2-pointcorrelation 
function, since the correlator is a 2-pointstatistic, by construc- 
tion, and so it becomes a measure of the underlying intrinsic 
2-pointcorrelations in the CMB (albeit via a specially selected 
subset of pairs of points on the matched circles). Therefore it is 
necessary to take into account possible variations in the underly- 
ing 2-pointcorrelation function with varying angular separation, 
which if not properly accounted for in simulations may lead to 
under(over)-estimation of the confidence level thresholds. 

Giv en that the concord ance best fit LCDM cosmological 
model ('Sper gel et ani2007h yields a very poor fit to the CMB 
data at large angular scales, due to lack of correlations in the 2- 
pointcorrelation function of the data with respect to the LCDM 
model at scales > 60°, and that the correlation statistic is sensi- 
tive to the details of the intrinsic 2-pointCMB correlations (and 
in particular to any large scale anomalies), we do not assume the 
LCDM model to help create our simulations. Instead, we take 
a model independent approach. As the CMB reference power 
spectrum in our GRF simulations of the expected signal, we use 
the reconstru cted power s pectrum from the three-year WMAP 
data (IHinshaw et al.l2007l) R Furthermore, we neglect the effects 
of cosmic variance, and only randomize the phases (and noise 
realizations) in our simulations. We remove the Cf=o,i (i-S- the 
monopole and dipole) components from our simulations. 
We use the same set of a;.,„s representing the CMB signal for 
a single simulation of the two datasets, followed by convolution 
with instrumental beam profiles. 

For each differential assembly (DA), we simulate the noise ac- 
cording its properties and scanning strategy (number of observa- 
tions per pixel in map) using uncorrected, Gaussian noise. 

The simulations are preprocessed in exactly the same way as 
the observational data. We neglect the impact of the (resolved 



In the high I end (noise dominated range) of the reconstructed 
power spectrum, the unphysical negative values are zeroed to have a 
zero contribution to the total variance of the map. This approximation 
has a negligible effect due to small statistical weight of the large 1 mul- 
tipoles, and large exponential Gaussian smoothing that we apply to the 
data. In practice, this approximation has a negligible effect on the vari- 
ance of the resulting simulation. Moreover, it can at most only make our 
analysis more conservative. 
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or unresolved) point sources which is negligible, since we ap- 
ply relatively large smoothing and use the Kp2 sky mask for the 
analysis. 

In Appendix A, we discuss the sensitivity of our results to 
the degree of smoothing, the sky mask applied, and the assumed 
statistical approach in greater detail. 



3. Statistics 

We describe our correlator statistics, parameter space, search op- 
timization and approach for assessing the statistical significance. 



3.1. Matched circles test 



As in ICornish et alj ( |2004|) and lRoukema et alj ( |2004|) we use a 
correlation statistic of the form 



S 



{Tfniimj) + {Tjiriirrij) 



(1) 



where the index i defines a set of all points in the "first" set of six 
circles related to the orientation of a fundamental dodecahedron; 
index j is the set of corresponding points along the matched six 
circles; and mi,mj are cut sky weights of the Kp2 sky mask, 
which can have a value of either for a masked pixel or 1 for an 
unmasked pixel. Clearly, perfectly matched circles would yield 
S = 1, which, due to non-zero noise contributions, is not possi- 
ble in reality. 

The dispersion of the correlation coefficient as defined in 
Eq. [T] is statistically enhanced in the small circles regime, due 
to the joint effect of the reduced number of points probing the 
matched circles as compared to larger circles, the accidental cor- 
relations of large (w.rt. the smoothing scale) flat fluctuations that 
happen to have similar (or opposite) large scale trends, as well 
as due to the fact that the rm.s. values necessarily shrink (down 
to zero in case of zero mean fluctuations) for circles of size 
comparable or smaller than the smoothing length. As shown in 
Sect, m this reduces the ability to robustly determine the degree 
of consistency or inconsistency of data with simulations, due to 
the finite accuracy of the S values and significant steepening of 
confidence-level contours in this regime. 

We note that the S statistic value would tend to unity, re- 
gardless of the shape of the underlying CMB fluctuations, as the 
monopole increases in the CMB maps. One could expect a simi- 
lar effect for the dipole component, for small circles. This effect 
would affect the simulations and the data to the same extent. The 
sensitivity of the test would however be significantly weakened 
and as such we remove the monopole and dipole components 
from the datasets for the analysis, and defer study of the impact 
of other small £ multipoles to appendix lAl 

3.2. Parameter space 

We perform a resolution-limited, full parameter-space search 
over the orientation of the fundamental dodecahedron, and over 
a limited range of the identified circles sizes of up to 20°. The 
parameters are defined as follows: l,b - galactic longitude and 
latitude of the first circle, g - the angle of rotation of the do- 
decahedron about the axis determined by {I, b), a - the angular 
radius of the matched circle, and s - the twist parameter defining 
the relative phase offset of the matched circles. 



We use the following parameter space: 

I G [0°,72°) 

b G [26.57°, 90°) 

g G [0°,72°) 

a G [1°,20°] 

s G {-36°,0°,36°} 



(2) 



The boundaries in {I, b) conservatively cover a larger region 
than the one twelfth of the sphere from which a "first" circle 
centre can be chosen non-redundantly. The range of angle g is 
72°, to cover all possible orientations of the fundamental dodec- 
ahedron for a chosen "first" circle centre. Values larger than 72° 
would yield the same set of 12 circle centres as a rotation by that 
angle modulo 72°. The interval in circle size a is chosen to be 
ro ughly symmetric and ce ntered about the 11° value suggested 
by Roukema et al] (20041) . The three twists s are chosen as in 
iRoukema et al.l (|2004|)" 

For all datasets we use the same resolution of 1° in probing 
the parameter space, except for the data with smoothing length 
A = 0.5°, in which case we use a resolution of 0.5°. 



3.3. Accuracy and search optimization 

The resolution of the data that we analyze is spatially constant 
and is limited by the finite pixel size, so circles of different sizes 
are probed by different numbers of pixels. 
As the parameter space of the search is large, it is important to 
consider the trade-off between the accuracy of the estimates of S 
(directly related to the number of pixels probing the underlying 
fluctuations) and the numerical computational time needed to 
obtain better accuracy. However, the speed of the search can be 
substantially increased, since the effective resolution of the data 
in our case is not limited by the pixel size, but rather by the 
smoothing length. 

In this section, we focus on the density of points (probing the 
fluctuations along circles in the sky) required to obtain a given 
accuracy in estimating S, and its dependence on the angular ra- 
dius of the circles a, a circle sampling density parameter r, the 
map resolution parameter ris, and the smoothing length proper- 
ties of the data. 

Maps smoothed with larger smoothing lengths have fewer 
significant, high-spatial frequency Fourier modes, and there is no 
need for fine sampling in order to fully encode the information 
content along the circles. Assessing the same level of precision 
for smaller circles also requires a smaller number of pixels than 
for larger circles. 

We perform a series of tests to determine the sampling den- 
sity required to achieve our desired accuracy level. The tests rely 
on measuring the speed of convergence to the "ideal" fiducial 
S^'^ value, derived using far more points in the circle than the 
number of available pixels along the circle in our datasetfl as a 
function of the increasing sampling density. 

We empirically model the circle sampling-density function 
in such a way that for a given r value parameter, and for a given 
smoothing length of the data, the accuracy of the resulting S 
values (i.e. the statistical size of the departure from the fiducial 
value) is approximately the same for all circle sizes (Fig. |2] left 
panel). We use the following fitted function: 



N, 



pix 



(3.40a[deg] + 76.85) (^) ( ^) (3) 



V256/ 



^ For all directions pointing inside a single pixel, the same tempera- 
ture value of that pixel is used. 
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Fig 2. Left paneL Convergence of 5* values to the "ideal" fiducial value S^'^ = S{r = 1000) as a function of resolution parameter 
r (a sampling density resolution parameter, defining the number of pixels to be used, to probe the CMB fluctuations along circles 
through Eq. ^) and as a function of circle size a. The A^pix(a) function shape for a given smoothing length is fitted linearly, so 
that the accuracy in S values was approximately constant for all considered circles radii. The assumed working precision level of 
AS" 0.01 is marked with thick horizontal line. For clarity, only AS* relation derived for data smoothed with Gaussian A = 0.5° 
is shown. Similar relations are obtained for the remaining three smoothing lengths. The average value (black thick line) from all 
tested circles radii is used to define the required value of r parameter for the circle search with data smoothed to 0.5°, in order to 
achieve the targeted accuracy on S value. Right panel: Average S convergence relations derived for data with different smoothing 
lengths A G {0.5°, 1°, 2°, 4°}, along with la error bars from 20 simulations. The intersections of these with the horizontal (black 
thick) line give the required values of r for each smoothing length in order to obtain the assumed working precision of AS" = 0.01. 



where A^pj^ is the number of pixels used for calculation of 5*, for 
a circle of angular size a, and for a map of resolution ?is. The 
resolution parameter r controls the sampling density. In prac- 
tice, we choose the closest, even integer as an Apj^ value for the 
calculations. This empirically-devised formula yields approxi- 
mately the same accuracy of derived values of S for all circle 
sizes (Fig.|2]left panel), and holds for all smoothing lengths. The 
aim is to find a value of r, for each smoothing length, which will 
provide sufficient accuracy. 

We therefore calculate the deviation AS{r) 



AS{r)^{\S,{r)~Sf\r^lOOO)\) 



(4) 



where the () averaging is performed over all curves derived from 
Eq.[3]for circle radii a[deg] G {3, 5, 10, 20, 30, 40}. 

We assume the working accuracy for S values to be AS* = 
0.01 throughout the analysis. This defines the required val- 
ues of the sampling density parameter r (Fig. |2] right panel) 
and the corresponding number of pixels to be used (Eq. O 
to achieve the targeted accuracy. For the smoothing lengths 
A[(ie(7] € {0.5, 1, 2, 4}, the required resolution parameter val- 
ues are r G {26, 18, 12, 8}. We use these values throughout the 
analysis with both the data and the simulations. 

3.4. Statistical significance 

In this section we discuss our statistical approach for quantifying 
the confidence intervals. 

Since our simulations simulate CMB fluctuations in an 
isotropic, simply-connected Universe, we test the consistency of 
the WMAP data with the null hypothesis that the CMB is an 
arbitrary realization of the GRF in a simply-connected space. 
We quantify the degree of consistency via S correlator values 
obtained from the data, and compared with those of simulated 
distributions from N^irn = 100 GRF simulations (Sect. |2]i. As 



an alternate hypothesis we choose the PDS topological model. 
The inconsistency of the data with the simulations, at high sig- 
nificance level, would then be considered as consistency in favor 
of the alternative hypothesis (PDS model). 

Since we are interested only in the highest positive S cor- 
relations, we build probability distribution functions (PDFs) 
of S^^^{a), the maximal value of the correlation S{a) found 
in the matched circle search in the parameter space {l,b,g,s) 
(Eq. |2]i, using N^i^^ = 100 simulation^ We probe the un- 
derlying PDFs of S^^^{a) at 8 different values of a, i.e. for 
a G {1, 2, 5, 8, 11, 14, 17, 20} in degrees. 

We reconstruct the confidence intervals [c(a), d{a)], for the 
68% and 95% confidence levels defined by the (cumulative) 

probability P of finding a GRF, simulated, CMB realization that 

crnax 
'data 



yields > 



PiSZ"" > ^dT)(«) 



c(a) 



(5) 



1 - 



E 



1=1, s 



max^^max 
sim,i — data 



1/^si 



where c{a) = min(5™^'')(a) and /(S""'''', a) is the MC probed 
PDF of the S""^'' values. 

We interpolate confidence interval contours for the remain- 
ing a values of the parameter space using 4th order polynomial 
fit. 



"* Statistically there are some small differences in the S values re- 
sulting from probing slightly different angular separations (arising due 
to different separations of pairs of points, for the same pair of matched 
circles, when calculated with two different phase twists: s — and 
— ±36), due to the dependence of the S value on the underlying CMB 
2-point correlation function. For such small twists, this is found to be of 
the same magnitude as the statistical error on the AS'(~ 0.01) for all 
considered circle radii. 
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In the next section we apply this procedure to the considered 
WMAP datasets and simulations and present our results. 

4. Results 

In Fig.[3]we present results of the all-parameter-space search for 
the WMAP first year ILC map (left panel), and the three year 
WMAP INC map (right panel). 

The signal at - 11° in Fig. 4 of iRoukema et al.l (|2004|) is 
reproduced and plotted with red crosses in Fig. [3] (middle-left). 

Clearly, it is not necessary to process large number of simu- 
lations to resolve high confidence level contours, since all the 
datasets are consistent with the simply connected space GRF 
simulations at a confidence level as low as about 68% at all 
smoothing scales. 

Is is easily seen that as the circle size shrinks to zero (a < 
2°), it is difficult to estimate precisely the significance of the 
detections since the CL contours steepen, while the accuracy of 
the S value determination is fixed at AS ^ 0.01. This effect is 
most severe for large smoothing scales, as expected. 

We note that the correlations S tend to increase in relation to 

the smoothing length applied to th e data. 

In particular, the signal reported in lRoukema et alj ( |2004|) is sen- 
sitive to increases in the smoothing length. While at the smooth- 
ing length of A = 1° there is practically no excess maximum in 
S-max^a 12°) for s = -36° relative to the 5™^'' values for 
s = 0° and s — +36°, on the other hand, the excess is clearly 
seen at the smoothing length of A = 4°, where its significance 
increases almost up to the 95% CL. 

The results for the three year INC data are consistent with the 
first year data, in the sense that no statistically important excess 
correlations are found. 

In addition to analysing our two primary datasets, we also 
carried out the following complementary searches. 

We completed an all-parameter-space search using the three 
year WMAP ILC data. We find that the excess correlation cor- 
responding to the hypothesized PDS model is weakened for all 
smoothing lengths (for s G {0°, +36°}), and basically indistin- 
guishable from the noise of what would be false positive detec- 
tions if we were to define the 68% confidence level as a detec- 
tion threshold. For s = —36°, we plot the S^^^{a) values for 
the three yeai" ILC data with a black line in the left column of 
Fig.H 

Our other complementary test was that we performed a 0.5° 
resolution all parameter space search, using first and three year 
INC data and did not find extra strong localized correlations. 
However since the computation time increases with the power 
of the increased resolution (i.e. increasing the resolution by a 
factor of two increases the calculation time by a factor of 2" 
where n = 4 is the number of parameters in parameter space) 
we haven't performed the significance analysis with simulations, 
and therefore we do not present these results. 

5. Discussion 

The analysis of the correlations derived from the data and pre- 
sented in the previous section finds no statistically-significant 
detections. The cross-correlations of the S values, obtained for 
different angular radii of the matched circles, were however ne- 
glected. It is of course faster to compute confidence intervals for 
a sparse parameter space and interpolate in between. However, 
the significance of any detections found this way (i.e. condition- 
ally to the a priori assumed circle radii) might be overestimated. 



compared to the case when all possible correlations were ac- 
counted for in the full covariance matrix analysis. In the present 
work, since we do not find any significant deviations from the 
null hypothesis (i.e. we do not find any strong outliers in the S 
correlations) in any individually probed value of the "a" parame- 
ter, we find no need for any further extensions to the significance 
analysis already pursued. 

We note that these cross-correlations are present not just in 
the data, but also in the Monte-Carlo simulations, so they affect 
the analyses of simulations and data to the same extent. 

Finally, we note that our statistical approach of considering 
only the maximal 5* correlation values could be altered to con- 
sider the full distributions of S correlations. However, we are 
especially interested in viable candidates for non-trivial topol- 
ogy (especially in the proposed correlation signal around angu- 
lar circle radii of 11°) and as such, the models with the largest S 
values are the best candidates. Given that there can be only one 
correct orientation of the fundamental dodecahedron and hence 
only one S correlation value corresponding to it (most likely the 
largest locally found S value), in the alternative way involving 
the full distributions of S correlations, the test would be heavily 
dominated by numerous values that will not be associated with 
the true topological correlation signal. As a result, the test would 
mostly measure a degree of consistency between the simulations 
and data with respect to underlying two-point correlations via 
the circles on the sky, rather than candidates for non-trivial topol- 
ogy. We therefore build and rely on the statistics of specifically 
selected (in the full parameter space search) S values, one for 
a given simulation at a given circle radius, to build our statis- 
tics and reconstruct the confidence thresholds. Although we are 
aware that relying on the distributions of maximal values of ran- 
dom variates may lead to asymmetrical distributions with en- 
hanced tails, we note that in the case of the S statistics, the 
possible values are by definition restricted to within the range 
[—1,1]. We also note that it is unnecessary to resolve high con- 
fidence level contours, since the data are consistent with simula- 
tions mostly to within ~ Icr confidence contour. 

6. Conclusions 

In iRoukema et al.l (l2004h it has been suggested that the shape 
of the space might be consistent with the Poincare dodeca- 
hedral space (PDS) model (Fig. This suggestion was due 
to an excess positive correlation in the matched circles test 
(ICornishetalj|l998bl) of the first year WMAP ILC map, how- 
ever the statistical significance of this excess was not specified. 

We have revisited those results and found consistent correla- 
tion excess corresponding to the same orientation of the funda- 
mental dodecahedron using independent software. 

We extended and updated the matched circles search with 
the WMAP three year ILC data and the three year foreground 
reduced, inverse noise co-added map, and tested these at three 
different smoothing scales (1°, 2°, 4°). 

We performed an analysis of the statistical significance of 
the reported excess, based on realistic and very conservative MC 
GRF CMB simulations of the datasets. 

We find that under "matched circles" tests, both the first and 
three year WMAP data are consistent with the simply-connected 
topology hypothesis, for all smoothing scales, at a confidence 
level as low as 68%, apart from the first year ILC data smoothed 
to 4°, which are consistent at 95% CL. 

Acknowledgements. The authors would like to thank the anonymous referee for 
useful comments and suggestions. BL would Hke to thank Naoshi Sugiyama for 
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Fig 3. Results of the search in the parameter space (see Sect. I3.2l i for the highest S correlations in the first year WMAP ILC map 
(left column) and three year INC map (right column) smoothed to A = 1° (top), A = 2° (middle), A = 4° (bottom). For clarity only 
the highest 72 S{a) statistic values are plotted for each of the three considered phase shifts: —36°, 0°, 36° marked with red (+), 
green ( x ) and blue (*) respectively, and separated by 0.2° offset for better visualization and comparison. The red crosses correspond 
to the PDS model. The thick solid line in the left column show the S™^^ values for a search in the three-year ILC data with a phase 
shift of —36 °. The 68% and 95% c onfidence level contours from N^^^ = 100 simulations are over-plotted. Clearly we reproduce 
the results of lRoukema et al.l (l2004l) . Most of the points with the highest S values in the range of 10° < a < 12° closely correspond 
to the solution depicted in Fig. [l] It is easily seen that much higher correlation coefficients would have been required in order to 
significantly reject the null hypothesis that the Universe is simply connected in favor of the PDS model alternative hypothesis. 
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model will be statistically imprinted onto the S values evalu- 
ated from simulated maps (generated according to the assumed 
model). This will lead to biases in estimates of the confidence 
level thresholds. Given that the nature of the large scale (mostly 
quadrupole and octupole) anomalies of the WMAP data w.r.t. 
concordance LCDM model is unknown, we assumed a model 
independent approach for generating simulations as described in 
Sect. |3] In this section we show the impact of biasing of S val- 
ues due to large scale uncertainties in the assumed CMB power 
spectrum. 

In order to quantify the impact of the underlying CMB power 
spectrum on the S correlator values, we perform the following 
exercise. 

We create two INC simulations of the WMAP three year data, 
for two different power spectra, yet keeping exactly the same 
phase information in both cases. 

One simulation is made using exactly the power spectrum recon- 
structed from the WMAP observations as in Sect.|2](i.e. neglect- 
ing cosmic variance effects). 

The other simulation is made using a random realization of the 
best fit LCDM model (Hinshaw et al. 2003) but with the same 
phase information as in the first simulation. We Gaussian smooth 
them to common resolution with beam of FWHM= 1°. This 
guarantees that the differences in S values will only be due to 
different underlying power spectrfl More specifically, since the 
LCDM model yields a good fit to the WMAP data for large /'s, 
the discrepancies will be only due to the large scale anoma- 
lies. Using 20 random orientations of the fundamental dodec- 
ahedron, we then calculate the average of differences of the S 
correlator values between the two maps. That is, we calculate: 
(AS*) = {Slcdm - SwMAp)- Hence, the positive (AS*) val- 
ues show the excess correlations that one would additionally get 
if the LCDM model was assumed, as compared to the "exact" re- 
alization of the reconstructed WMAP power spectrum (and vice- 
versa for the negative values). 

In Fig. lA.ll the (AS*) relation is shown with red line. 
Clearly, the strongest additional correlations appear for the 
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Appendix A: Power spectrum, cut sky and 
smoothing iengthi dependence on S 
correiations. 

A.1. S dependence on the CMB power spectrum 

As mentioned in Sect.|2]the S correlation value depends not only 
on the particular alignment of CMB fluctuations, but also on the 
underlying CMB 2-pointcorrelation function (or its Legendre 
transform - power spectrum), as it is also a two point statis- 
tic. Therefore any discrepancies of the data from the assumed 



-0.1 I ' ' ' ' ' ' ' ' 1 
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Fig A.l. Average difference (A5) = {Slcdai — Swmap) cor- 
relator values, as a function of circle size, statistically probing 
different regions of the underlying power spectrum (red line). 
With green line also the difference relation is plotted for the ori- 
entation of the dodecahedron depicted in Fig.[r| 



' The noise component is negligible at the considered smoothing 
scales of 1° 
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smallest circles, i.e. for the largest angular separations in the 2- 
pointcorrelation function, and are as large as ~ 0.15. (These 
scales are additionally contaminated due to large smoothing 
scales as compared to the circle size). While for the circle sizes 
of about a ~ 10° the effect is small, it is obvious that at smaller 
circle sizes the confidence level contours obtained from the sim- 
ulations performed according to the LCDM model would be too 
conservative than it is needed. Therefore, in this regard, given 
that the LCDM model is currently widely accepted, we consider 
our analysis to be very conservative. This plot includes the de- 
pendence of (A5) for larger circle sizes than those which we 
analyzed. And it is easily seen that for scales of about 40° the 
confidence level intervals would be underestimated. 

We therefore conclude that the simply connected space hy- 
pothesis would be consistent with the data at yet even smaller 
CL, if the LCDM model were assumed for the generation of 
simulations. 

A.2. S dependence on galactic sky cut 

In this section we show the effects of different galactic sky cuts 
on the resulting amplitude of the correlation signal. For this pur- 
pose we work with the first year WMAP ILC map, smoothed to 
a resolution of 2° without mask. 

We compute the the correlation statistic s for the ori- 
entatio n of the dodecahedron found in iRoukema et alj 
(12004 and confirmed in this work (i.e. for the do- 
decahedron with the following face centers {l°,b°) = 
{(252, 65), (51, 51), (144, 38), (207, 10), (271, 3), (332, 25)} 
and their opposites). We then apply different sky masks and 
show the S{a) relations in Fig. IA.2I We used three different 
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FigA.2. Effects of different sky masks on the S{a) correlation 
statistic. See text for mask abbreviations definition (Sect. IA.2I ). 
One sigma CL contour derived for the case of Kp2 is also plot- 
ted. 



is lower (most likely due to some residual foreground contami- 
nations of the Galactic plane). 

We do not show the galactic sky cut dependence for the INC 
dataset. This is because there are stronger foregrounds present 
than for the ILC dataset, so that it is not straightforward to pres- 
mooth the map without contaminating the "clean" regions of the 
sky by the Gaussian tails of smoothing kernel. 

A.3. S dependence on smoothing lengtli 

In this section we present the impact of different smooth- 
ing lengths used during map smoothing process on the S 
correlation values. We show the dependence using the dodec- 
ahedron orientation corresponding to the highest correlation 
value (S^^^ value) at a = 12° in Fig. [3] (middle-right) 
which corresponds to the dodecahedron with faces {l°,b°) = 
{(49, 51), (91, -2), (144, 40), (256, 64), (333, 23), (28, -10)} 
and their opposites. This closely matches the dodecahedron 
depicted in Fig. [T] We use the WMAP three year INC map and 
the Kp2 sky mask. 

In Fig. IA.3I we show the S{a) values for three different 
smoothing lengths used in the analysis. The characteristic trend 
towards higher correlation coefficients as smoothing length in- 
creases is clearly apparent, but the the 68% CL contours increase 
by about the same amount. 
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FigA.3. Effect of different smoothing lengths on the S{a) cor- 
relation statistic. We plot S{a) dependence for the three consec- 
utive smoothing lengths used in the analysis. Corresponding la 
confidence thresholds for each smoothing scale are also shown. 



sky masks: which we call "bgc", " Kp2", and "kp03". The first 
of these is the one that was used in iRoukema et al.l ( |2004|) (i.e. 
the galactic plane is masked for |6| < 2°, and points within 20° 
from the Galactic Center are masked). The second and third sky 
masks correspond to the Kp2 and the third year KpO sky masks 
dBennett et al.ll2003h . 

Clearly, the fine details on the sky masking in case of ILC 
maps do not have large impact on the resulting S correlation 
value, except for the "no mask" case where the correlation peak 



